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Abstract 



The detection of an afterglow, following a 7-ray burst (GRB), can 
be described reasonably well by synchrotron emission from a relativis- 
I tic spherical expanding blast wave, driven by an expanding fireball. 

' We perform detailed calculations considering the emission from the 

0^ ■ whole region behind the shock front. We use the Blandford McKee 

(—| . self similar adiabatic solution to describe the fluid behind the shock. 

Qh| Using this detailed model, we derive expressions for the peak flux, 

Q ■ and the peak frequency at a given observed time. These expressions 

provide important numerical corrections to previous, more simplified 
^ I models. We calculate the observed light curve and spectra for several 

magnetic field models. We show that both the light curve and the 
- spectra are flat near the peak. This rules out the interpretation of the 

r> I optical peak of GRB970508 as the peak of the light curve, predicted 

' by the existing fireball models. We calculate the observed image of 

a GRB afterglow. The observed image is bright near the edge and 
dimmer at the center, thus creating a ring. The contrast between the 
edge and the center is larger at high frequencies and the width of the 
ring is smaller. 



1 Introduction 

The detection of delayed X-ray, optical and radio emission, "afterglow" , fol- 
lowing a GRB is reasonably described by emission from a spherical relativistic 
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shell, decelerating upon collision with an ambient medium (Waxman 1997a, 
Meszaros & Rees 1997, Katz & Piran 1997, Sari, Piran & Narayan 1998). 
A relativistic blast wave is formed and expands through the surrounding 
medium, heating the matter in it's wake. The observed afterglow is believed 
to be due to synchrotron emission of relativistic electrons from the heated 
matter. The surrounding medium will be referred to as interstellar medium 
(ISM), though this may not necessarily be the case. 

At any given time, a detector receives photons which were emitted at 
different times in the observer frame, at different distances behind the shock 
front and at different angles from the line of sight (LOS) to the center of the 
GRB. The properties of the matter are different at each of these points, and 
so are the emissivity and the frequency of the emitted radiation. Early calcu- 
lations have considered emission from a single representative point (Meszaros 
& Rees 1997, Waxman 1997a, Sari, Piran & Narayan 1998). Later works in- 
cluded more detailed calculations. Synchrotron emission was considered from 
the shock front (Sari 1998, Panaitescu & Meszaros 1998), and monochromatic 
emission was considered from a uniform shell (Waxman 1997b). 

We consider an adiabatic hydrodynamical evolution, and use the adia- 
batic self similar solution found by Blandford and McKee (1976) for a highly 
relativistic blast wave expanding into a uniform cold medium, which we will 
refer to as the BM solution. We neglect scattering, self absorption and elec- 
tron cooling. Self absorption becomes important at frequencies much smaller 
than the peak frequency, and for slow cooling, electron cooling becomes im- 
portant at frequencies much higher than the peak frequency, so this should 
yield a good approximation for the observed flux around the peak. An anal- 
ysis of the spectrum over a wider range of frequencies was made by Sari, 
Piran and Narayan (1998). 

In section 2 we derive the basic formula for the observed flux from a 
system moving relativistically. In section 3 we consider synchrotron emission 
from a power law electron distribution, and calculate the observed light curve 
and spectra. We show that both the light curve and the spectra are flat near 
the peak. This causes difficulty in explaining the shape of the optical peak 
of GRB970508 (Sokolov et. al. 1997). We obtain expressions for the peak 
frequency at a given observed time, and for the peak flux. In section 4 
we consider three alternative magnetic field models. We show that the light 
curve and the spectra remain fiat near the peak in all the cases we considered. 

In section 5 we calculate the observed light curve and spectra from a 
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locally monochromatic emission. We consider a uniform shell approximation, 
which was calculated by Waxman (1997c), and compare the light curve and 
spectra to those obtained for the BM solution. We show that a uniform 
shell approximation yields results which are considerably different from those 
obtained for a more realistic hydrodynamics, and very different from those 
obtained for a realistic emission together with a realistic hydrodynamics. 

In section 6 we calculate the surface brightness, thus obtaining the ob- 
served image of a GRB afterglow. As indicated in previous works (Waxman 
1997b), Sari 1998, Panaitescu & Meszaros 1998), we obtain from detailed 
calculations, that the image appears brighter near the edge and dimmer near 
the center, creating a ring near the outer edge. At a given observed time, 
the contrast between the edge and the center of the image is larger and the 
width of the ring is smaller at high frequencies, while at low frequencies the 
contrast is smaller and the width of the ring is larger. 

2 The Formalism 

We consider a system that is moving relativistically while emitting radiation. 
We obtain a formula for the flux that is measured by a distant detector (i.e. 
D ^ L where D is the distance to the detector, and L is the size of the 
area emitting radiation). We use a spherical coordinate system and place 
the origin within the emitting region (i.e. at the source), while the z-axis 
points towards the detector (see Figure The detector is at rest in this 
frame, and so is the ambient ISM. We refer to this frame as the observer 
frame. Consider a small volume element dV = r'^drd^dcj) (where /i = cos^) 
and let H, be the energy per unit time per unit volume per unit frequency per 
unit solid angle, emitted by the matter within this volume in it's local frame 
(note that generally depends on the direction VL as well as on the frequency, 
place and time). We denote quantities measured in the local rest frame of the 
matter with a prime, while quantities without a prime are measured in the 
observer frame. Note that jv/^"^ is Lorentz invariant (Rybici & Lightman 
1979) and v' = 2/7(1 — Pfiv), where 7 and /5c are the Lorentz factor and 
the velocity of the matter emitting the radiation, respectively, and fiy is the 
cosine of the angle between the direction of the velocity of the matter and 
the direction to the detector, in the observer frame). The contribution to I^, 
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from this volume element is given by: 



dl = 3jy dr 

^"^V 7^(l-/5/i.)2^ 

(see Figure p. The contribution to the flux at the detector is dF^, = I^dfl 
where Q is the solid angle seen from the detector, and I^, includes all the 
contributions from different volume elements along the trajectory arriving at 
the detector from the direction Q simultaneously and at the time for which 
Fi^ is calculated. Consider a photon emitted at time t and place r in the 
observer frame. It will reach the detector at a time T given by; 

T = t-^, (2) 

where T = was chosen as the time of arrival at the detector of a photon 
emitted at the origin at t = 0. Using a = ry/1 — fi'^/D, we obtain: 



where j' = v' ■, ?, t), and later P' = P'{v' ., r, t). j' is taken in the direction 
Q'^ at which a photon should be emitted in order to reach the detector, and 
7, P, fiy should be taken at the time t implied by equation 

For a spherical expanding system, which emits isotropically in it's local 
rest frame, one obtains fiy = and Jl = PI^/Att, so that: 

Note that because of relativistic effects, a jet with an opening angle 6 > l/'j 
around the LOS can be considered locally as spherical (Piran 1994). 

In order to learn whether the radial integration is important, we calculate 
the observed flux from emission only along the LOS. We do this by consider- 
ing a situation in which at each point the photons are emitted only radially: 
Jl, = Pl^,S{Q' — Q'{f)). Note that the correct limit is obtained when the delta 
function in the direction of the emission is taken in the local frame. Since 
rfll = 7^(1 - (3^YdVt' (Rybici & Lightman 1979) we obtain that: 

6{n' - n\r)) = 7^(1 - (3^i)H{n - n{r)) . (5) 
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Figure 1: The contribution of a volume element dV to the flux observed by 
a distant observer is dFi, ^ I^dfl, where = j^ds. Since the observer is far 
away, the direction of emission in the observer frame is almost parallel to the 
z-axis. 

Substituting this / in equation ^ we obtain: 

Equation ^jis quite general, and includes integration over all space. In the 
case of GRB afterglow, radiation is emitted only from the region behind the 
shock front. The spacial integration should therefore be taken over a finite 
volume, confined by the surface of the shock front. We would therefore like 
to obtain an explicit expression for the radius of the shock R as a. function 
of /i = cos 6 for a given arrival time T. In the case of a shell moving with a 
constant velocity /5c, one obtains from equation ||: 

If one considers a constant arrival time T, this equation describes an ellip- 
soid, which confines the volume constituting the locus of points from which 
photons reach the detector simultaneously (Rees 1966). In GRBs, most of 
the matter is concentrated in a thin shell which decelerates upon collision 
with the ambient medium. When the deceleration of the shell is accounted 
for, the ellipsoid is distorted. The details of this distortion depend on the 
evolution of the shock radius R{t) (Sari 1998, Panaitescu & Meszaros 1998). 



5 



In this paper we consider an adiabatic ultra-relativistic hydrodynamic solu- 
tion, which imphes T oc i?"^/^, where T is the Lorentz factor of the shock. 
For this case, equation ^ yields: 

^=l-/i+V(8P) • 

The shape of the volume constituting the locus of points from which photons 
reach the detector simultaneously resembles an elongated egg (see Figure 0), 
and will be simply referred to as the "egg" . The side facing the observer (the 
right side in the Figures) will be referred to as the front of the "egg", and 
the side closer to the origin (the left side in the figures) will be referred to as 
the back of the "egg" . 

From hydrodynamic considerations, we expect the typical size of the re- 
gion emitting radiation behind the shock, to scale as: A cx R/T^ in the 
observer frame. Despite this fact, it is still important to consider the emis- 
sion from the whole volume of the "egg" , whose limits are given by equation 
^. To illustrate this we give a simple example. Consider a photon emitted 
on the LOS at a distance d = behind the shock (point A in Figure |^) 

at a time to in the observer frame, where 7 = T/\/2 is the Lorentz factor of 
the matter just behind the shock. This Photon will catch up with the shock 
front at a later time ti in the observer frame (at point B in Figure 0), and 
will arrive at the detector together with a photon emitted at that point at 
ti- From equation II we find that on the LOS: R = 167^cT, and so we obtain 
that R{ti)/R{to) = 17^/^ = 2. This shows that the emission comes from a 
substantial part of the volume of the "egg", and not just from a thin layer 
near it's surface. This is illustrated by the shaded region in figures ^ and |ll], 
which corresponds to a shell of width A = R/Aj"^ in the observer frame. 

At a given observed time T, the emission should be considered from the 
volume of the "egg" whose surface is described by equation |^. Taking this 
into account, it is simpler to calculate the flux at a given observed time T, 
using new variables y, x that depend on T (see Figures ^ and |^): 



R ^ 2 ^R-r 



y^l ' X-l + 167|(^), (9) 



where R = R{t) is the radius of the shock front, 7/ is the Lorentz factor of 
the matter just behind the shock and Ri is the radius of the point on the 
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the shock front 
at time t„ 




Figure 2: We see here the egg-shaped region from which photons reach an 
observer at a given time T. In all the figures of the "egg" the observer is 
located far to the right, and the symmetry axis is the LOS to the center 
of the GRB. In order to reach the observer simultaneously, photons emitted 
at different locations should be emitted at different times in the observer 
frame, according to equation 0. Photons that are emitted simultaneously 
in the observer frame along the dashed-dotted line ac, reach the observer 
simultaneously. Therefore ac represents an equal y contour line. The location 
of the shock front at this time of emission is indicated by the solid line abc, 
and it's radius R appears in the definitions of y, and x- A photon emitted 
at point A, at a distance of = behind the shock front and at a time to 
(both in the observer frame) caches up with the shock front at point B, at a 
later time ti in the observer frame. 
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Figure 3: Equal-y contour lines {y = 0.1, 0.15, 0.95) - vertical lines, 
and equal-x contour lines {Logio{x) = 0, 0.15, 0.3, 3) - curved lines. The 
horizontal line is the LOS to the observer, which is located far to the right. 

shock front, on the LOS from which a photon reaches the detector at a time 
T (see Figure 0). Since we expect the typical size of the emitting region 
behind the shock, to scale as: A oc the choice of x is natural to this 

problem. The exact form of x was chosen to suit the BM solution, discussed 
below. Using equation ^ we can express r, ^ in terms of y, x- 

r^Riy , (10) 

167/1/ 

where 7/ is the Lorentz factor of the matter just behind the shock, on the 
LOS. 

We would like to express equation ^ in terms of y, x- This will enable us 
to calculate the flux for the BM solution. This solution describes an adiabatic 
highly relativistic blast wave expanding into an ambient uniform and cold 
medium (Blandford & McKee 1976). In terms of y and x the BM solution is 
given by: 

n' = 4jfn,x''^^ , 7 = 7/X"'/' , e' = 4n,m,c^^'fX-''^'\ (H) 

where n' and e' are the number density and the energy density in the local 
frame, respectively, 7 is the Lorentz factor of the bulk motion of the matter 
behind the shock, rup is the mass of a proton, ni is the number density of 
the unshocked ambient ISM in it's local rest frame and 7/ = ■jiy~^^'^. For the 
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BM solution, one obtains (Sari 1997): 

7. = 3.65E,f nr^/^T-^// , = 5.53 x 10^' EUW^^TZcm , (12) 

where E^2 is the total energy of the shell in units of lO^^erg, ni is the number 
density of the ISM in units of cm~^ and T^ays is the observed time in days. 

For any spherically symmetric self similar solution we can define g{x) by: 
7^ = '~fjg{x), where g{x) describes how 7 varies with the radial profile. For 
the BM solution g{x) = X~^j ^-^id for a uniform shell g{x) = 1. Using the 
definitions above, we obtain from equation ^ after the change of variables: 

F(u T) = ^ Td r"'' d y''P'iT,u^{l-Pfx),y,x) .^3. 
Ji ^Jo y [I + y^8g{x)-' - X)]' 9{X) ' 

This formula for the flux takes into consideration the contribution from the 
whole volume behind the shock front, and will be referred to as the general 
formula. 

We would like to single out the angular integration and the radial inte- 
gration, in order to find out the different effects each integration has on the 
observed flux. 

In order to single out the effect of the angular integration we consider a 
thin shell of thickness A in the observer frame and take the limit A — > 0. 
Because of kinematical spreading we expect A to scale as: A cx -R/7^. Ac- 
cording to the definition of x, such a shell corresponds to a constant interval 
in X 5 i-6. the shell lies within the interval x ^ [l^Xmax] for some Xmax- The 
limit A ^ corresponds to taking a delta function in x • ^{x " !)• 

In order to single out the effect of the radial integration, we changed 
variables in equation ^ from r to x- Fo^' the BM solution we obtain: 



3 Synchrotron Emission 

According to the fireball model, a highly relativistic shell moves outward and 
is decelerated upon collision with the ambient ISM. This creates a relativistic 
blast wave that expands through the ISM and heats up the matter that passes 
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through it. The relativistic electrons of the heated material emit synchrotron 
radiation in the presence of a magnetic field. 

We now consider the synchrotron emission at a certain point, in which 
the values of n', 7 and e' are given by the hydrodynamic solution. In order to 
estimate the local emissivity, we assume that the energy of the electrons and 
of the magnetic field at each point are a fixed fraction of the total internal 
energy at that point: e'^i = e^e' , = e^e'. We assume that the shock 
produces a power law electron distribution: N{'je) = K^e~^ for 7e > ^min 
(for the energy of the electrons to remain finite we must have p > 2). In the 
Figures for which a definite numerical value of p is needed, we use p = 2.5. 
The constants K and 'jmin in the electron distribution can be calculated from 
the number density and energy density: 

7mm = (- ^ I 2 ' K^iP- ^Wl^in ■ (15) 

\p — 1 J n'meC^ 

Assuming an isotropic velocity distribution, the total emitted power of a 
single average electron (i.e. with 7e = (7e)) is given by: 

Kav-t^TCP^^efe'^ , ile) ^ (16) 

6 nrrieC^ 

(Rybici & Lightman 1979), where ctt is the Thomson cross section and e'^ = 
S^/Stt, where B is the magnetic field (in the local frame). Although we refer 
to the magnetic field in the local frame, throughout the paper, we make an 
exception and write it without a prime. 

The synchrotron emission function (power per unit frequency) of a sin- 
gle electron is characterized by P^^e oc u^^^ for frequencies much smaller 
than the electrons synchrotron frequency, and it drops exponentially at large 
frequencies. The typical synchrotron frequency, averaged over an isotropic 
distribution of electron velocities, is given by: 

where nie and are the mass and the electric charge of the electron, respec- 
tively. We approximate the emission of a single electron as P^^e oc v^l^ up to 
the electrons typical synchrotron frequency, where we place a cutoff in the 
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Figure 4: The local emissivity from a power law distribution of electrons, 
emitting synchrotron radiation. The solid curve represents exact synchrotron 
emission with an isotropic electron velocity distribution, while the dashed 
curve represents the approximation we used for the local emissivity. The 
differences in the local emissivity tend to smear out when integration is per- 
formed over the whole volume behind the shock front. 

emitted power. We normalized the emission function so that the total power 
emitted by a single electron equals that of an exact synchrotron emission. 

Under these assumptions, we obtain after integration over the power law 
electron distribution, that the spectral power per unit volume (in the local 
frame) at any given point is: 



where u'^^^ = i^'synilmin) is the synchrotron frequency of an electron with 
the minimal Lorenz factor at that point. Since the emitted power at each 
point peaks at z^^j^, this frequency can be looked upon as the typical emit- 
ted frequency around which the emitted power is concentrated. Although 
this emission function was obtained by approximating the spectral emission 
of each electron as having the shape of the low frequency tail (oc z/^/'^), the 
spectral power for the whole electron distribution resembles that obtained 
for an exact synchrotron emission and an isotropic electron velocity distri- 
bution. The solid curve in Figure ^ represents the local emissivity from an 




(18) 
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exact synchrotron emission of an isotropic distribution of electrons (Rybici & 
Lightman 1979, Wijers & Galama 1998), while the dashed curve represents 
equation ITB| with: 



(3p- 1) Z/^j^„((7e)) 

where a factor of 0.88 was added to improve the fit to the exact synchrotron 
emission. In our calculations we use the local emission represented by the 
dashed line in Figure |. This local emission differs from that of an exact 
synchrotron emission near the peak, by up to ~ 45%, and is only slightly 
different above or below the peak. Note that differences in the local emis- 
sivity tend to get smeared out, when the contribution to the observed flux 
is integrated over the whole volume behind the shock front. We expect that 
considering exact synchrotron emission should somewhat lower the peak flux 
and the peak frequency, and make the light curve and spectra more rounded 
and fiat near the peak. We evaluate that the peak flux would be lower by 
about 30%. Since there are only slight differences in high and low frequen- 
cies, there should hardly be any effect on the value of the point where the 
extrapolations of the power laws at high and low frequencies meet (see Table 

!)• 

In this section we use the BM hydrodynamical solution. The combination 
of the realistic emission function described above, and a realistic hydrody- 
namical solution should yield a reasonable approximation for the observed 
flux. 

The results are presented as a function of the dimensionless similarity 
variable = v/vt-, where vt is defined as the observed synchrotron frequency 
of an electron with 7e = 7mm just behind the shock on the LOS: 

UT ^ vmm{y = X = 1) = 1.54 X lO^^v/m (^^j ef e^^Esf r^-^Z/Hz 

(20) 

where z is the cosmological red shift of the GRB. Our results can therefore 
be looked upon, with the proper scaling of the logarithmic x-axis, either as 
the spectra at a given observed time or as the light curve at a fixed observed 
frequency. 

Similarly, we express the observed fiux in terms of a "standard" fiux Fq 
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Figure 5: The observed flux from synchrotron emission of a power law 
electron distribution for the BM solution. The different curves stand for 
emission only along the LOS, emission from an inflnitely thin shell, and 
emission from the whole volume behind the shock. The emission from the 
whole volume is very similar to that obtained for an inflnitely thin shell, due 
to a coincidence that rises from our choice for the magnetic fleld model (see 
section |^). 



deflned by: 



where P^^^^{Ri) and v'synini) are the total power and synchrotron frequency 
of an average electron at |/ = x = 1; respectively, Pu,peak is an estimate of 
the peak spectral power of an average electron, and N[Ri) = AnRfni/S is 
the total number of electrons behind the shock at the time t in the observer 
frame for which R{t) = Ri. Allowing for cosmological corrections, Fq is given 
by: 

Fo = 8.24 X 10^(1 + z)e]i'E,2ny' ' ^^^^ 

where di is the luminosity distance of the GRB. 

The results for the observed flux are presented in Figure |. The flux 
arriving from the LOS (the dashed curve) peaks at (ppeak — 0.98, and is only 
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slightly rounded near the peak. The flux arriving from an infinitely thin shell 
(the dotted curve) peaks at (ppeak — 1-92 and is quite rounded and flat near 
the peak. The flux arriving from the whole volume behind the shock front 
(the solid curve) peaks at 4>peak — 1-88 and is flat and rounded near the peak, 
quite resembling the flux from an infinitely thin shell. 

Our best prediction for the flux measured by a distant detector is ob- 
tained when we assume a realistic hydrodynamical solution (the BM solu- 
tion), a realistic radiation emission and use the general formula (see the solid 
curve in Figure |]). We would now like to examine it more closely. The 
curve looks quite flat near the peak, and we attempt to demonstrate this 
feature in a quantitative manner. If one compares the peak flux obtained 
by extrapolation of the power laws, obtained at low and high frequencies, 
to the "actual" peak flux, one obtains that it is larger by a factor of 1.53: 
Fextr = 1 -53^1, max- In ordcr to further estimate the flatness of the curve 
near the peak, we define 0+ and 0_ by Fj,(0+) = F^{(j)J) = 0.5F^^rnax, where 
(^+ > ippeak > (p-- We obtain that 0_ ~ O.10peafc and 0+ ^ lO^peafc, so that 
~ 100. In a similar manner, we define T+ and T_ for a given observed 
frequency, by F„{T+) = Fi,{T^) = 0.5F,,,max, where T+ > Tpeak > T_. One 
obtains: T^/T^ = (0+/0_)^/'^ ~ 20 (note that T_|. and T_ are frequency 
independent). 

The frequency at which the observed flux peaks at a given observed time 
T is given by: 



The main dependence on p lies in the factor (p — 2)^/(p — 1)^, but there is 
also a weak dependence of (ppeak on p. The value given here for (ppeak is for 
p = 2.5, and it decreases with increasing p (it is higher by 24% for p = 2, 
and lower by 12% for p = 3). The maximal flux is given by: 



Here there is a much weaker dependence on p. The value stated here is for 
p = 2.5, and it is smaller by 11% for p = 2 and larger by 6% for p = 3. 

These results can be best understood by looking at the "egg" that con- 
stitutes the locus of points from which photons reach the detector simul- 
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J^peak = (ppeakJ^T = 2.9 X lO^Vl + 



^T4Eli%-SR. . (23) 



F,, 




(24) 
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Figure 6: Equal observed synchrotron frequency contour lines, at a given 
observed time, for the BM solution. Each line represents a different value of 
0. The bold contour line is = 1. The lines are separated by a constant 
logarithmic interval in 0: ALogiQ{(f)) = 0.05 (the inner contour lines, which 
cross the LOS, represent < 1). The shaded region represents a finite shell 
of thickness A = in the observer frame. 

taneously, and mapping upon it the typical emitted frequency in the ob- 
server frame. Figure |^ depicts the lines of equal I'syn- From the elec- 
tron distribution we obtain that (7e) = 'jminip ~ ~ 2) and therefore: 
^min = ijP ~ 2)^/(p — lYi'syn- For this rcasou, the frequency contour lines 
that are depicted in Figure ^ also represent equal Vmin lines. 

There are two opposing factors that determine the shape of the observed 
synchrotron frequency contour lines, that are depicted in Figure The 
shift from the local frame to the observer frame, v' = — jSu), causes 
the frequency to decrease as one moves away from the LOS. However at 
earlier emission times, and at locations closer to the shock front, the typical 
synchrotron frequencies in the local frame are higher, and so is the Lorentz 
factor of the matter. This tends to increase the observed frequencies of 
photons that were emitted earlier (i.e. from the back of the "egg") and 
closer to the shock (i.e. closer to the surface of the "egg"). 

The result of these two opposing effects, for the BM solution, is that as 
one goes backwards along the LOS, to earlier emission times, the observed 
synchrotron frequency for a constant observed time is almost independent of 
X- UsyniLOS) oc x~^^'^^- This explains the result for the flux arriving from 
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the LOS (the dashed curve in Figure ||), namely that the peak flux is obtained 
at a frequency just slightly lower than vt'- 4>peak — 0.98, and the light curve is 
only slightly rounded near the peak. The fact that < along the LOS 
accounts for the fact that the peak flux is obtained at (ppeak < 1- The fact 
that the decrease in z/^m as one goes back along the LOS is very moderate, 
means that in order for the emitted radiation to be concentrated around a 
frequency substantially lower than ut, one has to go very deep in the radial 
profile along the LOS to r <^ Ri. This implies that one gets far from the 
shock front (to x ^ 1), and therefore the contribution obtained to the total 
flux is small. For this reason (ppeak is very close to 1. 

The observed flux from an infinitely thin shell peaks at (ppeak — 1.92, 
and the light curve (or spectra) is quite rounded and fiat near the peak (see 
the dotted line in Figure ||). This result can be understood by following 
the surface of the "egg" (see Figure ^. z/^m increases as one goes to the 
back of the "egg" (i.e. to earlier emission times) along it's surface, and it 
does so much faster than it decreases when one goes back along the LOS. 
Therefore one gets a substantial contribution to the flux at > 1, before one 
gets too far back in the shell, where the total contribution to the flux drops 
considerably. This explains why the peak flux for an infinitely thin shell is 
obtained at a frequency significantly higher than ut, whereas for the LOS it 
is obtained at a frequency only slightly lower than z/y. 

The contribution to the flux from an infinitely thin shell, as a function 
of the radius of the shell, for different values of is shown in Figure 0. The 
different curves can be looked upon as representing either different observed 
frequencies at the same observed time, or the same observed frequency at 
different observed times. One can see that for (p < 1 the contribution to the 
flux peaks a.t y = 1 (the LOS), while for > 1 it peaks at smaller values of y. 
One can obtain the average radius contributing to a given (p in the following 
way: 



lody- ay 



(y) = Z : ^ • (25) 



{y){(f)) is shown in Figure |. For (p < 1 one obtains that (y) is constant: 
(jj) = 0.82, while for > 1 there is a decrease in (y) until for ^ 1 it reaches 
the limiting value: (y) = 0.65. The same infinitely thin shell approximation 
was considered by Panaitescu & Meszaros (1998) and by Sari (1998). The 
values quoted by Panaitescu & Meszaros for (y) in the two limits of high 
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Figure 7: The contribution to the observed flux as a function of y = R/ Ri, 
for an infinitely thin shell. The different curves stand for LogiQ^cj)) = 
—1.5, —1.25, 1.75 and = ^peaifc ~ 1-92. For < 1 the contribution 
peaks at the LOS {y — 1), and for > 1 the contribution peaks at a lower 
value of y as increases. 
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Figure 8: The average normalized radius {y) contributing to the observed 
flux at a given 0, for an infinitely thin shell, {ij) ^ 0.82 for < 1, and it 
decreases as increases above 0=1, until it reaches the limiting value of 
{y) ~ 0.65 for large 0, at around 4> ~ 100. 

and low frequencies are probably erroneous, and they are not consistent with 
the value they quote for the bolometric {y) (this inconsistency occurs only 
for the case we consider in this paper, of an adiabatic evolution of the shell, 
and neglecting electron cooling). The values obtained by Sari are slightly 
different than ours, since he considered emission from a shell consisting of a 
given element of matter. 

4 Alternative Magnetic Field models 

Although the hydrodynamic profile is given by the BM solution, the structure 
and profile of the magnetic field are less clear. So far we have assumed that 
everywhere the magnetic field energy is a fixed fraction of the internal energy: 



Since not much is known about the origin or spacial dependence of the mag- 
netic field, we consider now two alternative models for the magnetic field. 
This can serve, in a way, as a test for the generality of our results, and 




(26) 
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it will expose any hidden fine tuning caused due to our previous magnetic 
field model, if one existed. We assume that each matter element acquires a 
magnetic field according to equation ^ when it passes the shock. The two 
alternative models are obtained by assuming that the magnetic field is either 
radial or tangential, and evolves according to the "frozen field" approxima- 
tion. 

Consider a small matter element, which passes the shock at a time to in 
the observer frame. Just after it passes the shock it possesses a magnetic 
field -Bo (in the local frame) given by equation |2^. We consider a cubic 
volume Vq = (in the local frame) with one face perpendicular to the radial 
direction. According to the BM solution, at a later time t this matter element 
will be at X = (Rii) / Rito)y , and it will occupy a box of a size Lrad = X^^^L 
in the radial direction, and a size = x^^^L in the two tangential directions 
(0 and 0). One also obtains that Bq = Bjx^^^, where Bj is the magnetic 
field at the front of the radial profile, just behind the shock, at the time 
t. We consider two possibilities for the direction of the magnetic field at 
to: a radial and a tangential magnetic field. Brad and B±. Our previous 
"equipartition" model will be denoted simply by B. In both new cases the 
"frozen field" approximation implies that the magnetic field will remain in 
the same direction, while it's magnitude changes in the following way: 

Brad = BfX-'^' , B = BfX-''^'' , B^ = BfX-'. (27) 

Using the general formula (equation |^), we obtain now the observed flux 
for Brad and B^ (see Figure P). Some of the features are summarized in 
Table 1. Since the emission from the shock front (x = 1) is identical in all 
three cases, the results obtained for an infinitely thin shell in the previous 
section, including the dotted curve in Figure |, are still valid here. This 
implies that the difference between the various curves in Figure |^ arises from 
the radial integration. Our previous model for the magnetic field led to an 
almost negligible effect of the radial profile (see Figure Now we find that 
generally, the effect of the radial integration is comparable to that of the 
angular integration, as the flux for B^ad and -B^ are substantially different 
than for an infinitely thin shell. 

From equation we can see that Brad implies a larger magnetic field 
than B, and therefore a larger total emission and higher emission frequencies, 
while B±^ implies a smaller magnetic field than 5, and therefore a smaller 
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Figure 9: The observed flux for tfircc different magnetic field models. Bj-ad = 
B = B± on the shock front, while behind the shock front Brad > B > B±. 



TABLE 1 



Features of the Light Curve and Spectra for 
Different Models for the Magnetic Field. 



Magnetic 
Field Model 


(ppeak 




0- 


0+ 


0+ 
4- 


T+ 


Fo 


4^extr 


B± 


1.67 


0.178 


0.157 


15.9 


101 


21.7 


0.273 


4.12 


B 


1.88 


0.205 


0.187 


18 


99 


21.4 


0.313 


4.66 


Brad 


2.86 


0.284 


0.241 


29 


120 


24.3 


0.454 


7.30 



Table 1: 0peafc, -Pl/.maa:, -^0, ct, 0-, and T+ are all defined in the text. 

The last two columns refer to the point where the extrapolations of the power 
laws at high and low meet. 
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total emission and lower emission frequencies. This is consistent with the 

results in Table 1, namely: 4>peak,rad > (ppeak > (ppeak,± and Fi,^rnax,rad > 

Fu,max > F^,max,±- An important feature that appears in all the magnetic 
field models we considered, is the flatness of the peak: ~ 100 — 120 

and T+/T_ ~ 21 - 24 (see Table 1). 



5 A Uniform Shell Approximation 

In this section we consider a locally monochromatic emission. We consider a 
case where all the electrons at a certain point have the same Lorenz factor, 
which we chose to be the average Lorenz factor of the electrons at that 
point, given by equation We take the angle a between the velocity of 
the electrons and the direction of the magnetic field in the local frame, to be 
the average value obtained for an isotropic velocity distribution: sin a = 7r/4. 
Under these assumptions, the emission from this point is at a single frequency, 
given by equation |T^. The results are still expressed by the similarity variable 
(f) = u/ut, where in this section ut can be obtained from equation ^ by 
dropping the term involving p. The emitted power per unit volume per unit 
frequency in the local frame is given by: 

Pi, = n'P^^^Jiu' - z.;j , (28) 



where is given by equation |l| . 



We now turn to calculate the observed flux due to this emission within 
a thin shell of matter, of finite thickness A in the observer frame. This case 
was calculated numerically by Waxman (1997c). He took: A = (R/j'^, which 
corresponds in our notation to: XmaxiC) = 1 + 16C- We use equation p!3|, and 



integrate only up to XmaxiC)- Substituting g{x) = 1 in equation |I3|, We 
obtained an analytic solution, and present it in Figure We compare this 
both with a locally monochromatic emission from a BM profile (see Figure 
1^), and with synchrotron emission from a BM profile, obtained in section ^ 
The flux obtained from a finite shell with a "cut" BM profile, is very 
similar to the flux from a full BM profile, shown in Figure |1^, except that 
there is a cutoff at < (1 — 160^^/^"^, due to the finiteness of the shell. The 



place of this cutoff is shown in Figure |T2|, for a few values of C, including 
( ^ 0. Figure |10] depicts the observed flux from a uniform shell, for the 
same values of (, and is considerably different. This large difference can be 
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Figure 10: The observed flux from a uniform shell, for different values of 
the shell thickness in the observer frame: A = (R/'y'^. The value of at the 
peak is a function of the shell thickness (see equation |29| ). 



understood by looking at Figures ^ and |TT], which depict the contour lines 
of equal observed emission frequency (which will be referred to as frequency 
contour lines) for the BM solution and for a uniform shell, respectively. 

We now turn to Figure |12], which depicts the flux from a BM profile with a 
locally monochromatic emission. Normalizing according to the maximal flux, 
the flux from the LOS and the flux from an infinitely thin shell, practically 
coincide with the flux calculated from the general formula, for < 1 and 
> 1, respectively, and therefore are not shown separately. The flux peaks 
at (ppeak = 1 and drops sharply for < 1, while exhibiting a more moderate 
decrease for > 1. The frequency contour lines for the BM solution are 
nearly parallel to the LOS and their length within the thin shell (or up to 
a constant value of x) drops quickly to zero for < 1 (see Figure H). This 
explains the quick drop (and eventual cutoff, for a finite shell) of the flux for 
< 1, and the fact that the peak flux is located at = L If one follows 
the shock front (i.e. the surface of the "egg"), one finds that the observed 
frequency increases as one moves from the LOS to the back of the "egg" (i.e. 
to smaller radii). This explains why we find a cut-off in the flux arriving 
from an infinitely thin shell, at frequencies smaller than ut (0 < 1). Since 
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Figure 11: Equal observed synchrotron frequency contour lines, at a given 
observed time, for a uniform shell. The lines from right to left, represent: 
(f) — 5^/^, 5^/'^, 5^/^, .... The typical synchrotron frequency is larger than vt 
everywhere (everywhere > 1), and increases as one goes to the back (i.e. 
to the left) of the "egg" along the LOS, to smaller radii. The shaded region 
represents a finite shell of thickness A = i?/47^ in the observer frame. The 
bold frequency contour line that touches the back of the shell on the LOS 
at point C, represents the frequency for which the flux from a uniform shell, 
considered in this section, peaks. 
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Figure 12: The observed flux from a locally monochromatic emission at the 
typical synchrotron frequency, using the BM solution. Since oc z/T^/^, this 
figure represents (with the correct scaling) both the light curve at a given 
frequency, and the spectrum at a given observed time. The peak flux is 
obtained at = 1, and the flux drops sharply for smaller 0. If we consider 
only a slice of a BM profile, this introduces a cutoff at low 0. This cutoff 
occurs at points A through D, for ^ — 0, C = 1/64, ( = 1/16 and ( = 1/4, 
respectively, where the shell thickness in the observer frame is A = (R/'y'^. 



there is a sharp drop in the contribution to the flux as one goes back in the 
radial profile (away from the shock front), the flux at > 1 where there 
is a contribution from the front of the radial profile, is similar to the one 
obtained for an infinitely thin shell. For < l,on the other hand, there is no 
contribution from the very front of the radial profile (near the shock front), 
and therefore the contribution from back the radial profile becomes apparent, 
and we obtain a result similar to the one obtained for the LOS. 

The observed flux for a uniform shell is shown in figure |10[ The shape 
of the light curve and the spectra, the shape and sharpness of the peak and 
the value of (ppeak, all change considerably for different values of the shell 
thickness (. The big difference compared to the results for a BM profile, 
can be understood by looking at Figure |TT], which depicts the frequency 
contour lines for a uniform shell. The frequency contour lines are nearly 
perpendicular to the LOS and higher frequencies come from the back of the 
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■'egg" (i.e. from smaller radii). For this reason the flux vanishes for < 1 



(see Figure |T0]). For every value of ( there exists a critical frequency 0c(C) 
for which the frequency contour line touches the back of the shell, exactly on 
the LOS. This is demonstrated by the bold frequency contour line in Figure 



TT| , which represents 0c for ^ = 1/4, and touches the back of the shell exactly 



on the LOS at point C. For < 0c(C) there is a sharp drop in the length of 
the frequency contour lines within the thin shell, that causes a sharp drop 
in the flux. For > 0c(C) there is a moderate change in the length of the 
frequency contour lines within the thin shell, and the flux drops (though more 
moderately) for > 0c(C) as well. This causes the peak flux to be exactly 
at 0c(C)- From our analytical solution for a uniform thin shell we obtained 
a simple analytical expression for the time (or frequency) of the peak flux as 
a function of 

^peakiC) = MO = (1 + 160'/' , (29) 

The resulting spectra and light curve for a uniform shell are qualitatively 
different from those obtained when the hydrodynamical radial profile is con- 
sidered. For a uniform shell, the rather arbitrary choice of ( determines the 
time of the peak flux for a given observed frequency, and the peak is sub- 
stantially sharper, with a smaller width at half maximum. The difference is 
even bigger if the results for a uniform shell are compared to those obtained 
for a realistic emission and a realistic hydrodynamic solution (the solid curve 
in Figure ^ . 



6 The Observed Image 

We turn now to calculate the observed image, as seen by a distant observer. 
We consider the BM hydrodynamic solution. Substituting g{x) = X^^ i^i 
equation |13| we obtain: 

8Rf PV\ ^ 
dF. = -^^l^J^.dydx . (30) 

We calculate the surface brightness (energy per unit time per unit frequency 
per unit area perpendicular to the LOS) at a given arrival time T. The 
distance of a point from the LOS is given by: 

R^=rsme^ Riy^l - /x^ ^ ^^^y-^y^ . (31) 
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^ = 0.1 ^ = ^max=''-88 ^ = 10 




Figure 13: The observed image of a GRB afterglow, for different values of (j). 
The different images can represent either the image at different frequencies 
at the same observed time, or the image at the same frequency at different 
observed times. In the latter case, one should also remember that the size 
of the observed image increases with the observed time as R±^max oc T^/^. A 
larger value of represents a higher observed frequency or a later observed 
time. The image is brighter near the edge and dimmer near the center, and 
the contrast is higher at large values of 0. 



The maximal value of R± is obtained on the surface of the "egg" (where 
X = 1), implying y = 5~^/^. The observed image at a given observed time 
T is restricted to a disk of radius R±,max — 4:.14:'jicT around the LOS. Using 
equation |12|, we obtain an explicit expression for the BM solution: 

R±,max = 3.91 X lO'^^Eli'n^'/'T'.ilcm . (32) 

We calculate the surface brightness within this disk, and find it useful to 
work with the variable: x = R±/R_L^max- The differential of the area on this 
disk is given by: 

dS± = 2TTR±dR± = 2'kR\ ^^^xdx . (33) 



Using equation we obtain from equation ^ 



y 1 F'y^x , 



dS^ 7r\Dj cT{l + 7xy'^y 

from which we obtain after integration, the surface brightness as a function 
of R±. In this expression, x should be taken as a function of y for a given x, 
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TABLE 2 



Features of the Observed Image. 



Magnetic 
Field 
Model 


SB(x = 0)/SB^,, 


< 1 


0> 1 


<p < 1 


0> 1 


•^peak 


Ax 


■^peak 


Aa; 




0.32 


0.034 


0.95 


0.263 


0.96 


0.151 


B 


0.34 


0.039 


0.94 


0.296 


0.95 


0.178 


Brad 


0.39 


0.065 


0.93 


0.403 


0.90 


0.290 



Table 2: . We define Ax = x+ — a;_, where x+ > a;_ are the values of 
X = R±/ R±,max at which the surface brightness (SB) drops to half of it's 
maximal value. Ax is an estimate for the width of the bright ring that 
appears on the outer edge of the image. Xpeak is the value of x for which the 
maximal SB is obtained, and it's values indicate that the SB peaks near the 
outer edge of the image. It is evident that the contrast between the center 
and the edge of the image is considerably larger for high frequencies (or late 
times) than for low frequencies (or early times). 



according to equation The limits of the integration over y are determined 
from the condition x > 1- 

We calculated the surface brightness for several values of 0, and as before, 
the different images can be viewed as either the observed images at a given 
observed time at different observed frequencies, or as the observed images at 
a given frequency at different observed times. The observed images are quite 
similar for the different magnetic field models we considered, and therefore 
we present detailed results only for B (figures 13 and H), and summarize 
the features of the images obtained for Brad, B and B± in Table 2. 

The observed image is bright near the edge and dimmer at the center. At 
high values of (high frequencies or late observed times) the surface bright- 
ness at the center is only a few percent of the maximal surface brightness, 
which is obtained near the edge (x ~ 0.93 — 0.95), and a bright ring is clearly 
seen with a sharp edge on the outer side and a more gradual decrease in the 
surface brightness towards the center. At < 1 (low frequencies or early 
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Figure 14: The observed surface brightness as a function oi x = R±/ R±,max, 
for different vahies of 0. The surface brightness peaks near the outer edge 
of the image, drops sharply towards the edge (a; = 1) and decreases more 
gradually towards the center. The normalization preserves the correct ratio 
between the absolute surface brightness, for the different values of 0, which 
decreases both for low and for high values of (f). 



28 



observed times) the surface brightness at the center is around 0.32 — 0.39 of 
the maximal surface brightness, and though the image is brighter near the 
edge, the center should be visible as well (see the left image in Figure |T3[). 

The transition in the distribution of the relative surface brightness along 
the observed image, between the limiting cases of small and large 0, occurs 
over one order of magnitude in 0: 1 < < 10. The relative surface brightness 
hardly changes beyond this region (though the absolute surface brightness 
drops both for low and for high values of 0, due to the drop in the observed 
flux). 

The fact that at a given observed time, the relative surface brightness 
at the center of the image is smaller for high frequencies than it is for low 
frequencies, may be intuitively understood: for high frequencies the emit- 
tance decays in time, and is therefore lower at large radii, which correspond 
to small values of x oc R±, which are located at the center of the image. 



7 Discussion 

We have derived a formula for the observed flux from a relativistically mov- 
ing system. Using this formula, we considered the emission from the whole 
volume behind a spherical relativistic blast wave, expanding into a cold uni- 
form medium. We have considered an adiabatic evolution and used the self 
similar hydrodynamical solution of Blandford & McKee (1976) to describe 
the matter behind the shock. 

We have considered synchrotron emission from a power law distribution 
of electrons, using a realistic local emissivity, which depends on the local 
hydrodynamic parameters. We have obtained expressions for the frequency 
at which the observed flux peaks at a given observed time, and for the value 



of the peak flux (equations ^ and 24). The value we obtained for the peak 



frequency at a given observed time is a factor of ~ 1.8 smaller than the value 
obtained by Sari, Piran & Narayan (1998), a factor of ~ 2.5 (for p = 2.5) 
smaller than the value obtained by Wijers & Galama (1998), and a factor 
of ~ 35 (for p = 2.5) smaller than the value obtained by Waxman (1997b). 
The large difference from the result obtained by Waxman is partly due to the 
fact that he took the peak frequency to be at the synchrotron frequency of 
an average electron i^syniile)), while the peak frequency is actually obtained 
much closer to I'synilmin), i-e. the synchrotron frequency of an electron with 
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the minimal Lorentz factor. The value we obtained for the peak flux is a 
factor of ~ 1.9 (for p = 2.5) smaller than the value obtained by Wijers & 
Galama (1998), a factor of ~ 5.1 larger than the value obtained by Waxman 
(1997b), and a factor of ~ 6.5 smaller than the value obtained by Sari, Piran 
& Narayan (1998). 

We have shown that both the light curve and the spectra arc flat near 
the peak. The flux at a given observed time drops to half it s maximal 
value at around one order of magnitude from the peak frequency, on either 
side. The flux at a given observed frequency drops to half it's maximal value 
at about a factor of ~ 5 before and after the time of the peak flux. This 
result was obtained for all the magnetic flcld models we considered, and it 
therefore seems to be of quite a general nature. Our calculations were made 
using an approximation for the local emissivity that is obtained from an exact 
synchrotron emission. We expect that without this approximation the values 
of the peak flux and the peak frequency should be somewhat lower, and the 
light curve and the spectra should be even more rounded and flat near the 
peak. We estimated that these effects should be small. 

In contrast to the flatness of the peak, discussed above, GRB970508 dis- 
played a sharp rise in the optical flux, immediately followed by a power law 
decay (Sokolov ct. al. 1997). Sari, Piran & Narayan (1998) considered 
both radiative and adiabatic evolution of the blast wave, and found that the 
steepest rise in the flux occurs before the peak, for an adiabatic evolution 
and slow cooling of the electrons, which is the case discussed in this paper. 
This steepest rise is T"-*^/^, and as we have shown in this paper, the rise in the 
flux decreases as the peak is approached, and the peak itself is quite flat. We 
obtained T+/T_ ~ 21 — 24, which indicates a flat peak, while GRB970508 
displayed T+/T_ < 3. This rules out the interpretation of the optical peak of 
GRB970508 as the peak of the light curve, predicted by the existing flreball 
models. It therefore appears that another explanation is needed. 

We have considered a locally monochromatic emission, at the typical 
synchrotron frequency, obtained from the relevant hydrodynamic parameters. 
We have shown that when this emission is applied to the BM solution, the 
observed flux peaks at the observed emission frequency of the matter just 
behind the shock, on the LOS {(ppeak — 1)- There is a sharp drop in the flux 
below this frequency, and a gradual decrease above it. When this emission 
is applied to a uniform shell (Waxman 1997b) the results change drastically. 
In particular the location of the peak flux depends critically on the width of 
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the shell ( (see equation |2^) , which is chosen quite arbitrarily. 



The image of a GRB afterglow looks like a ring, even when emission is con- 
sidered from the whole volume behind the shock front. Similar results were 
obtained for more simplified models (Waxman 1997b, Sari 1998, Panaitescu 
& Meszaros 1998). The image is bright near the edge and dimmer at the 
center. The contrast in the surface brightness between the center and the 
edge of the image is larger for high frequencies (optical and x-ray) and the 
ring is narrower, while for low frequencies (as long as self absorption is not 
significant) the contrast is smaller and the ring is wider. Radio wave-bands 
can be considered as "low frequencies" for the first few months. The best 
available resolution is obtained in radio frequencies, and the afterglow of a 
future nearby GRB {z ~ 0.2) might be resolved. This theory predicts that in 
early times, when Uradio < ^peak the image should appear as a relatively wide 
ring with a relatively small contrast, while for later times where I'radio > ^peak 
(as long as the relativistic regime is not exceeded) the image should appear 
as a narrow ring and posses a large contrast. 

The observed image has been calculated, considering emission only from 
the surface of the shock front (Sari 1998). This yielded a surface brightness 
diverging at R±,max- This divergence is an artifact of the assumption that 
the radiation is emitted from a two dimensional surface. Other features of 
the image, such as the difference between high and low frequencies, are quite 
similar in the more simplified analysis. 
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